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__  ABSTRACT 

This  paper  developed  out  of  an  attempt  to '"u" nd er S tci nd  refeuits  of  Gelfand- 
Dikii  [1]  and  P.  Moerbeke  (unpublished  version  of  [2])  in  a unified  way.  We 
provide  a natural  abstract  setting  for  understanding  the  symplectic  structure 
involved  in  both  results.  The  setting  is  an  orbit  in  the  dual  algebra  of  a 
group.  We  also  discuss  a natural  trace  functional  for  formal  asymptotic 
pseudo-differential  operators.  In  addition  we  discuss  so-called  Lenard 
recursion  relations  inherent  in  these  structures. 
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SIGNIFICANCE  AND  EXPLANATION 


It  is  one  thing  to  write  down  the  equations  governing  the  behavior  of  a 
physical  system  in  time;  it  is  quite  another  thing  to  solve  them.  Historically, 
an  important  step  forward  in  explaining  the  behavior  of  mechanical  systems  was 
made  when  it  was  realized  that,  for  some  systems,  certain  quantities  like  energy 
and  angular  momentum  were  independent  of  time,  i.e.  invariant.  This  is  a tremen- 
dous help  in  solving  the  equations  both  analytically  and  numerically,  since  the 
invariants  are  equivalent  to  first  integrals  of  the  equations  of  motion,  and 
constrain  the  class  of  functions  within  which  we  need  look  for  solutions. 

It  turns  out  that  in  many  problems  of  physical  interest  involving  partial 
as  well  as  ordinary  differential  equations,  it  is  possible  to  find  a denumer- 
able number  of  quantities  that  are  invariant  as  time  evolves.  This  paper  gives 
a unified  method  for  finding  invariants  for  a class  of  equations  that  includes 
crystal  lattices  (the  Toda  equations) , stochastic  processes  (Kac-Moerbeke) , and 
water  waves  (Boussinesq  and  Korteweg-deVries) . In  addition,  the  mechanical 
structure  of  these  systems,  as  discovered  by  respectively  Gelfand-Dikii , and 
P.  Moerbeke  are  given  geometrical  significance. 
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The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MR C,  and  not  with  the  author  of  this  report. 
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ON  THE  RELATIONSHIP  OF  SOME  RESULTS  OF  GELFAND-DIKI  I AU!>  u.  .".  I.M  i r'i  , 

AND  A NATURAL  TRACE  FUNCTIONAL  FOR  FORMAL  ASYMPTOTIC  PSEUDO-DIFFER! NTT AI  : ERA TOP . 

M.  Adler 

1 • Introduction 

This  paper  developed  out  of  an  attempt  to  understand  the  results  of  Gelfand-Dikii  [1] 
and  P.  Moerbeke  (unpublished  version  of  [21),  in  a unified  way.  The  above  authors 
respectively  discovered  generalizations  of  the  Korteweg-deVr ies  equation  and  Toda 
system  equations,  which  are  completely  integrable  Hamilton  systems  whose  equations  of 
motion  are  expressible  in  terms  of  the  Lax-isospectral  equation. 

The  crucial  observation  in  both  cases  is  that  the  relevant  symplectic  structure 

in  these  problems  is  the  orbit  symplectic  structure  of  Kostant-Kirillov  [12,18].  We 

* 

understand  that  for  the  Toda  system,  this  fact  was  known  to  Kostant  . We  observe  that 
this  symplectic  structure  plays  a role  in  the  n-dimensional  Euler  spinning  top  problem 
of  V.  Arnold  [11],  along  with  an  isospectral  equation,  as  was  discovered  by 
L.  Dikii  [15].  In  addition,  quotient  symplectic  structures  through  Hamiltonian  group 
actions  [12),  of  which  the  Kostant-Kirillov  construction  is  but  an  example,  have  been 
seen  to  play  a role  in  the  Moser-Calogero  systems,  see  [13,14]. 

In  the  case  of  the  generalized  Toda  systems,  the  relevant  group  G is  the  group 

of  lower  triangular  matrices  with  nonzero  diagonal  elements.  We  identify  it's  dual 

. ★ 

algebra  with  the  upper  triangular  matrices  through  the  trace  form,  and  then  the 

orbit  Hamiltonian  phase  space  0^  of  interest  is  of  the  form 

eft  = {[u’1AU]+|u  e G},  A e £*  , 

where  [B]  + denotes  the  matrix  formed  from  B by  setting  its  lower  trianqular  entries 
equal  to  zero.  We  shall  wish  to  place  certain  conditions  on  A,  see  Section  2. 

* 

Personal  communication  from  J.  Moser. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-75-C-0024  and  the 
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For  the  case  of  the  Generalized  Korteweq-deVr ies  equation,  the  relevant  group  G 

is  the  formal  pseudo-differential  symbols  of  negative  type,  translated  by  the  identity 

* 

element  1,  whose  dual  C we  may  identify,  throuqh  a trace  form  to  be  specified,  with 

the  differential  symbols  of  nonneqative  type , which  are  identified  wi  th  differential  operators . 

Then  the  Hamiltonian  orbit  space  of  interest  is  of  the  same  form  as  0,  above,  where 

A 

the  operation  ( )+  shall  denote  the  natural  projection  of  a pseudo-differential 

symbol  onto  the  nonneqative  symbols.  For  the  case  of  the  Korteweq-deVr ies  equation, 

2 

A = -p  + q,  an  orbit  9^  is,  rouqhly  speakinq,  specified  by  the  condition 
J qdx  - constant . We  also  point  out  that  althouqh  the  formula  "in  coordinates"  for  the 
correct  symplectic  structure  of  the  Generalized  Korteweq-deVr ies  equations  appears  in 
111,  it  comes  about  throuqh  strictly  computational  procedures,  and  it  is  our  purpose 
to  place  the  formula  in  its  proper  geometric  sett i ng . This  also  necessitates  redoing 
parts  of  (1]  in  a less  ad  hoc  procedure  than  was  employed  in  [1],  so  as  to  give  complete 
and  consistent  proofs,  and  moreover  simplify  proofs  in  general  that  appear  there. 

We  amplify  the  previous  paragraph.  The  algebra  of  pseudo-differential  operators 
has  a subalgebra  consisting  of  those  symbols  which  have  an  asymptotic  description  at 
One  may  think  of  taking  this  description  as  a set  of  normal  corrdinates  for  the  sub- 
algebra, since  one  identifies  symbols  modulo  their  behavior  on  compact  sets.  This 
suggests  that  one  should  abstract  the  algebraic  content  of  the  above  situation,  and  just 
work  in  a formal  setting.  Namely  we  work  with  formal  Laurant  series  in  a variable  £ 
over  some  formal  differential  ring,  with  the  multiplication  rule  inherited  from  pseudo- 
differential operator  theory.  For  motivation,  one  thinks  of  the  formal  Laurant  series 
as  being  in  a real  variable  £ embedded  in  a complex  neighborhood  of  00 , or  more 
precisely  as  a series  germ.  We  then  define  the  trace  functional  (which  is  of  course 
commutative)  as  the  equivalence  class  of  the  coefficient  of  the  £ ^ term,  where  we 
identify  two  elements  precisely  if  their  difference  is  a total  derivative  of  a ring 
element.  The  trace  functional  makes  possible  the  identification  of  a cotangent  bundle, 
and  from  then  on,  the  Korteweq-deVr ies  case  proceeds  in  perfect  analogy  with  the  Toda 
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a point  which  we  wish  to  stress.  The  method  also  works  automatically  fox 
the  se  ] f -a»i  ioi  nt  case,  which  m«*r<  Iv  mv-lv-s  a sliahtly  different  coord  1 n 1 zat  ion  of  the 
formal  Laurant  series  then  the  previous  case,  or  for  the  cast*  where  the  relevant  rinq  is 
composed  of  square  matrices  (see  [l*Jl)  whose  entries  are  themselves  differential  rinq  elements. 

In  the  latter  cast-,  we  must  take  as  the  trace  functional,  the  equivalence  class  of 
the  matrix  trace  of  the  coefficient  of  the  £ * term. 

We  point  out  that  if  the  semi-lower  triangular  matrices  have  a place  in  the 

Fourier  theory  of  say  Z,  or  Z , the  integers  and  the  integers  modulo  n,  respec- 

n 

tively.  There  place  is  analogous  to  the  role  of  the  formal  asymptotic  pseudo- 
differential operators  of  negative  type.  We  discuss  this  briefly  at  the  end  of  Section  2. 

We  conclude  in  Section  4 with  a generalized  set  of  so-called  Lenard  relations,  which 
provides,  at  least  in  theory,  an  alternate  method  of  computing  all  the  relevant  quantities 
in  the  theory,  and  in  addition,  in  the  case  of  self-adjoint  operators,  the  relations 
contain  formal  spectral  information  concerning  the  operators  studied.  These  relations 
are  essentially  a consequence  of  the  law  of  exponents  for  operators.  There  are  a 
number  of  conjectures  concerning  these  relations.  We  note  that  Bill  Symes  has  discussed 
recursion  relations  in  [8]/  and  the  author  in  [5]. 

I wish  to  thank  Bill  Symes  for  many  stimulating  discussions.  I am  also  indebted 
to  J.  Moser  who  provided  the  stimulus  and  encouragement  for  this  line  of  work,  along 
with  the  clarification  of  some  points.  Tn  addition,  my  thanks  extend  to  Joel  Robbin 
and  C.  Conley,  for  helpful  suggestions,  and  at  whose  seminar  at  the  University  of 
Wisconsin,  Madison,  these  results  were  first  presented.  We  also  mention  that  Bill  Symes 
has  a generalization  of  the  Gelfand-Dikii  theorem  in  the  setting  of  this  paper. 
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1 1 . The  Generalized  Tod a Systems 

As  motivation,  and  of  interest  in  its  own  right,  we  firsf  discuss  the  Toda  system 
and  its  generalizations,  which  display  the  relevant  structures  in  both  examples. 

To  this  end  we  introduce  the  group  G of  lower  triangular  matrices  with  strictly 
nonzero  diagonal  entries.  Its  Lie  algebra  £ is  just  the  lower  triangular  matrices, 

h 

and  we  identify  £ , the  dual  of  £,  with  the  upper  triangular  matrices  (see  l IB]) 

via  the  trace,  ( *>,  functional  and  associated  pairing  ( •,*>, 

(2.1)  < A>  = tr IA.  . ] = l A.  . , < E,F>  = < EF>  . 

ij  n 

★ * 

Since  g e G acts  on  £ via  conjugation,  it  naturally  acts  on  £ e £ through 

* 

duality,  and  with  our  given  representation  of  £ its  action  is  described  by 


(2.2) 


-1, 


g : £ ■*  [g  £ g]  + 


where  [ 1 denotes  the  projection  operation  of  setting  all  terms  below  the  diagonal 

equal  to  zero.  We  denote  an  orbit  of  this  action  by  8 = 6^*  = { [g  £ g]+|g  e G}. 

★ 

From  (2.2),  the  tangent  space  of  8 at  £ is  described  by 

(2.3)  T02* {£*)  = {[£,£*]  + 1 £ e £}  , 

and  the  natural  symplectic  2-form  u,  of  Kostant-Kir illov  [121  , associated  with  the 
orbit  space  8 is  defined  via 

(2.4)  <*>([£*,  £xl+,  (£*,£21  + ) <£*)  = <£*,[£1,£21>  = < [£* , ^1  +,  £.,  > . 

If  we  let  A = [A. .]  be  the  running  variable  on  £ , and  H = H(A)  = H([A  .1)  be  a 
ri  13 

* 

function  on  £ , we  may  uniquely  define  H , the  gradient  of  the  function  H = H(A) 
with  respect  to  < , >,  via 

(2.5) 


X (H)  = (X,H  >,  H f ( , 
A A 


where  X is  a vector  field  on  £ , and  hence  also  identif  able  as  an  element  of  £ . Then  since 

0)([A,-H  ] , [A,£  I ) = < [A.  J£  ] ,H  > = [A,£,l  (H)  , 

At  2+  2+A  2t 

we  have  that  Hamiltons  equations  with  H = H(A),  associated  with  the  symplectic 

structure  of  (2.4)  is,  from  the  definition  w(X  ,Y)  = Y(H), 

H 


(2.6) 


A " XH=  [HA'fl), 
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We  note  that  the  H e (h| [A,Hfl)t  = 0,  for  all  A f fl form  an  algebra,  namely 
the  algebra  which  characterize  the  invariants,  or  constants,  of  the  orbit 
0 = 0 . The  description  of  an  orbit,  9„,  consists  in  finding 

C 

these  constants.  We  also  note  from  (2.4,6),  that  the  Poisson  bracket  associated  with 
«,{.},  is 
(2.7) 


(H(1)  ,H<2)  }(A)  ; u(X  . ,X  ) (A)  = <A,[H(1),H(2>]>  . 

H 1 H<2)  A A 


We  now  specialize  these  considerations,  and  for  that  we  introduce  some  notation. 
Let  the  shift  operator  £ , acting  on  R°  be  defined  by 

(2.8)  (£v)  . = v , v = (v  , . . . ,v  )T  e Rn  , 

l l+l  1 n 


where  we  define  v.  = 0 if  i f {l,2,...,n}. 

l 

In  addition,  we  associate  with  a « Rn,  the  multiplication  operator  a*  , namely 


(2.9) 

With  these  definitions,  we  now  define 


(a  • v)_.  = ajvj»  3 = 1,2, ...,n  . 


/->  -.ns  A _ r„|„_  r „ (i)  r i (i)  , (i)  (i)  (n-l-i)  „ . (i)  , „n-, 

(2.10)  A . = 1 A [ A — l a • £ , a = (an  ,a,  , ,a;  ,0,0 0),a  e Rl 

k , l 


k<i<j 


0 1 


Note  a'J1  = A,  For  example,  if  A = J a^£^,  B = £ £ .^(b^*  )«  then 


,(j) 

i 


i,i+3 


i>0 


i>0 


<A,B>  = £ (a.,b.),  with  ( , ) the  scalar  dot  product  in  R , ( , defined  by 

i>0  n 

(2.1)  . In  the  future,  we  shall  omit  the  dot  in  a • £ when  there  is  no  possibility 

★ 

of  misunderstanding.  Clearly  £ - A . g C a „ , and  £ = 4 . Along  with  the 

-n,0  -n,0  0, n 

grading  inherent  in  the  specification  of  the  A . 's,  we  have  the  natural  projections 

K r 3 

P.  into  the  A 

k,u 

* 

We  shall  restrict  A e J.  , such  that  A e 4 ,0<m<n,  and  we  observe  from 

0 ,m  — 

(2.2)  that  B € 4_  implies  0„  C 4 , and  also  [H„,B]  e A , which  we  indicate 

0,m  * B 0,m  B -n,m 

by  replacing  the  subscript  + in  (2.6)  with  a superscript  m,  i.e. 


(2.11) 


A = (H  ,A]  . 

A 


Since  from  (2.11),  [H_,A]  only  depends  on  the  part  of  H„  contained  in  A »,  we 

A A -m, o 
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may  as  well  restrict  H to  only  be  a function  of  the  A . *s  such  that  0 < i - i < m 

i / J ~ 

Note  that  clearly  An  is  not  an  orbit,  for  if  D C 0 , trB  = trD.  This  is  a 
0 , m B 

consequence  of  H = trA,  = I,  implies  (HA»A]m  =0  in  (2.11),  and  so  trA  is  an 
orbit  invariant.  It  is  easy  to  see  that  in  qeneral  there  must  be  more  invariants  and 
we  shall  not  discuss  this  point,  but  refer  the  reader  to  f 1 8 1 , which  discusses  this 
problem  fully  for  the  special  case  of  orbits  of  maximum  possible  dimensionality. 

We  compute  (2.11)  in  coordinates,  using  the  notational  (2.8) - (2. 10) . We  have  that 


i 

L 


A - l a,/.  Hfl  = l C3  ~y 
k=0  j =0  3a  ^ 


, where 


hence  by  (2.11), 


A = l a(k,Sk  - [VA,m  = l l £ 

k=0  v=0  j>0 

and  so  Hamilton's  equations  are 

a(v>  = A = y A . -A 

t<t+V  j>ol3Vj,t  t-],v+t  t, 


(v+j) 

(V+j)  . [rv  3H  ) 

[-  L<j> 

• 

[ 1 

3H 

] 

A ^ , | J,  M A ✓ + ✓ M 

kt,t+v+j  3A 

J t+v, 

t+v+ j J 

, u ^ v ^ m,  u ^ t ^ n 

For  example,  if  we  set  m = 1,  A.  . = b.,  A.  . , = a.,  then  the  above  expression 

i.i  i 3,3+1  3 

yields  for  (2.11)  , 

b.  = (a.  ,H  - a.H  ),  a.  = a.(H  - H,  ) , 

l l-l  a.,  la.  r lb.  b., 

l-l  l l l+l 

hence  since  {f,h)  = X F,  we  have 
H 

{f,h)  = Z (a.  ,H  - a.H  )F.  + Z a.(H  - H,  )F 

l-l  a.,  l a . b.  lb.  b . , a . 

l-l  li  l l+l  l 

the  well  known  Poisson  bracket  of  the  Toda  system  [2) . In  this  case  it  is  only  necessary 

to  impose  the  condition  E b.  = constant,  to  specify  an  orbit,  assuming  none  of  the 

a^'s  are  zero.  The  property  of  a^  being  zero  is  an  orbit  invariant. 

We  now  make  an  important  observation  due  to  P.  Moerbeke  in  an  unpublished  version 

+ 0 - + 

of  12],  but  first  some  notation.  We  write  every  matrix  M = M + M + M , with  M 

the  strictly  upper  triangular  part,  M°  the  diagonal  part,  etc.,  which  brings  us  to: 

+ T 

Theorem  1.  (P.  Moerbeke) . Let  L be  the  real  symmetric  matrix  A + (A  ) . 

Then  if  H = Hj.(A)  = tr  f(L)  = <f(L)>,  Hamiltonian  equations,  A = [HA,A]m,  imply 


that  I.  satisfies  the  Lax  isospectral  flow  equation 

(2.13)  L.  = fP,L],  with  P = P = f(L)+  - f(L)~  . 

Moreover  this  implies  that  (H_,H  } = 0 for  ail  polynomial  f,  g,  i.e.  the  H-  a 

f g ^ f 

in  involution  with  respect  to  { , }. 

Proof . We  qive  a 'functional*  version  of  P.  Moerbeke's  proof,  the  unpublished 
version  of  [ 2 J . In  the  course  of  the  proof  we  shall  show  the  necessary  fact 
[P,L]  e A 

-m,  m 

We  compute,  using  the  notation  of  (2 .8) - (2 .10) , 


911 9 < f (L)  > 


>a(i)  = ,a(i>  = <f'(L>'  = 


<f ' (L)  , £ 1 + £-1>,  1 < i < m. 


<f' (L)  , £ > , i = 0 


Hence  if  f(L)=  [ (f  1 (L)  ) (l)£  1 + J £ '-1  (f 1 (L)  ) ( • ’ , ^IL—  = 2 (f  • (L) ) (i)  - 6 . „ (f ' (L))  (i)  , 

i=0  i=-n  »aU) 

i = 0,1, ...,m,  and  thus 

(2.14)  H = [ 2f ' (L)  ~ + f ' (L)  °] 

A -m 

where  [ ] denotes  the  natural  projection  into  A 

-m  -m , 0 

We  now  use  the  fundamental  observation  [f'(L),L]  = 0,  observing  0 = [f'(L),L]  = 
[f'(L)  + f ' ( L)  ^ + f'(L)  ,L|  j hence 

[P,L]  = [f'(L)+  - f'(L)~,L]  = [L,f'(L)°  + 2f'(L)~]  , 


[L,f’(L)°  + 2f'(L)"]  = [A  + l",  f ' ( L)  ° + 2f'(L)~]  , 

and  so  [P,L]  = [A,f'(L)°  + 2f'(L)“l  = (A,f'<L)°  + [2f'(U]  l"1,  (and  thus  by  the 

+ + -m 

skew-symmetry  of  P,  [P,L]  e A ),  and  by  (2.14))  = |A,H  ]m. 

-m , m A 

We  thus  have  shown 

m 

A = [A,Hfl)  implies  L = [P,L],  H = Hf  . 

Since  L = |P,L],  P = Pf  when  H = 

H = <g(L)  > = <g'  (L)  • L>  = <g'  (L)  , [P,L]  > = <P, [L,g 1 (L)  ] > = 0 , 


since  (L,g'(L)]  = 0,  but 


dt  g 


H = {H  , H_},  and  so  we  have  proven  the  statement  con- 
g (J  f 


cerning  the  involuti-  » of  the  H^'s.  This  concludes  the  proof  of  Theorem  1.  We  refer 


the  reader  to  theorems  analogous  to  Theorem  1,  for  other  finite  systems,  in  [9],  [10], 

There  is  more  information,  namely  recursion  relations,  or  so-called  Lenard  relations 


to  be  gleaned  from  the  companion  identity  to  the  fundamental  relation  [f'(L),L]  = 0, 

k 


i +k  i 

namely  the  stronger  statement,  L = L 


I/',  which  we  shall  discuss  in  Section  4, 
but  for  the  moment  we  shall  go  on  in  the  next  section  to  the  differential  operator  case. 
In  preparation  for  the  next  section  we  mention  immediate  generalizations  of  the  above 
discussion. 


Remark  1.  One  generalization  is  that  one  may  consider  the  a 


(i) 


in  A = l 

0<  i<m 


. (i) 


to  be  n-periodic,  i.e.  £ acts  on  R , where  v e R means  v = ( ...,v  ,,v  ,v  , ...  ) 

n n -101 


and  v is  an  n-periodic  vector,  i.e.  v 


(i) 


j+n 


/.  for  all  j,  and  a e R . As 
3 n 


before.  A,  L,  may  be  represented  by  an  n x n matrix,  and  all  the  proofs  given  above 
apply  to  this  case,  with  slight  modifications,  provided  we  don't  think  in  terms  of 


matrices,  but  operators,  for  instance  A = J a^'p,  where  we  require  that  m [t!  . 

i=0  2 


An  equally  interesting  generalization  is  when  £ acts  on  R , where  v e R means 

v = ( . . . ,v_^,Vg,v  , . . . , )T,  and  A = j a^£*,  a ^ e R . In  the  later  case  the 

i=0 

(formal)  Lie  group  and  Lie  algebra  are  respectively 


G = { £ a^i*£i|a^i'  e R , a.^°'  > 0 for  all  j},  £ = { £ 

i<0  -1  i<0 


(i).ii  (i)  », 

i £ [a  e R } , 


and  we  have  for  the  dual  of  £ , 


* - { I a 

i=0 


fOii  (i)  oo 

C |a  € r , N < oo} 


We  finally  motivate  the  constructions  of  this  section,  and  especially  the 
next  section,  by  a brief  discussion  of  Fourier  theory  on  finitely  generated  abelian 
groups.  For  a more  complete  discussion#  we  refer  the  reader  to  [16] . Since  such  a 


group  can  be  decomposed  as  a direct  sum  of  copies  of  Z,  Z , where  Z is  the  integers, 

n 


Z the  integers  {0, 1, . . . , n-l}  with  modulo  n addition,  it's  sufficient  to  discuss  Z,  Z 
n r 


1 
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We  first  discuss  Z , with  g f Z , i.e.  g e fo,  i , . . . , n-1  1 . The  character 
n n 

qroup  of  Z , Z are  the  functions  from  Z into  the  complex  numbers,  which 

n n n 

multiply  by  pointwise  multiplication,  having  the  properties 

xfg^  • x(q2)  = x(gx  + g2>.  Ixtgii  = 1 < 


note  X (-g)  = X (q)  t with  * being  complex  conjugation,  and  <(1)  = 1.  Clearly  the 

2Tti 

k 

, k e {0,1 , . . . ,n-l } , and  thus 


character  group  Z^  = (x  = x ^ I ( j ) 


3 n 

V ak = e 


Z^,Z^  are  identical  as  abstract  groups.  Note  we  may  change  our  viewpoint  and  think 

of  X,  v {j)  as  a character  of  Z . If  we  define  the  inner  product  on  functions  of  Z , 
( • ) n n 

r»  * 

= > f . (g)f - (g)  , then  clearly  (x  . ( • ) # X,,  ( ' ) ) = n6  , and  viewing  X,.><i> 

n g6  z 1 2 Ik  Ik 

n 

as  characters  in  Z we  have  by  the  above  duality  (x,  % (j),x,  , (k))--  = n6 . These 

n (*)(*)  Z ik 

n 

two  identities  immediately  yield  the  Fourier  decomposition  of  functions  on  Z^  and 

the  Plancheral  identity,  namely  if  f is  a function  on  Z , we  have 

n 

r - -1 

f = 2,  f(j)x->  f(j)  = n (f,x.)/  i.e.  f(*)  is  a function  on  Z , 

je  Z 3 3 n 

n 

l|f||  \ = <Cf)z  = I |f(i)  |2  = n||  f ||  | = n(f,f)£  =n  l |f(j)|2. 

n n ie  z n n je  Z 

n n 

If  we  identify  the  functions  f f F(Z  ) with  n periodic  vectors  f in  R , by 

n n 

the  recipe  (f) . = f(j),  then  we  may  define  the  shift  operator  £ : F(Z  ) h>  F(Z  ), 
j n n 

through  ( f ) ] . = f(j+l)].  Clearly  the  characters  of  z are  a complete  list  of  the 

3 n 

eigenfunctions  of  £ which  have  the  value  1 at  j = 0,  i.e.  C (Xj)  = aj*j'  from 

00 

which  we  conclude  if  we  decompose  an  operator  B acting  on  F(Z^)  ~ Rn  into 

B — [ b^  • C ^ = P(£),  where  • denotes  pointwise  multiplication,  that 
jfZn 


B(XW)  = I ^b(:))x  = { l aV3’}  ’ (X.)  = P(a  ) 


and  so 


je  Z ” " je  Z 

n n 


B(f)  - B<  l f(j)x  ) = l f(j)  • (I  aV)  • X 1 = l f (j)P(a .)  ■ 

je  z J je  z k 3 3 je  z 
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1 


hence  we  have 


(2.15) 


B<X,>  = P(a  ) ' xfc.  B(f)  = l f (j)  P (a .)  ' X • < if  B = P (£  ) 


The  above  situation  is  to  be  compared  with  the  Fourier  theory  of  S . Here 

(x)  I , n e z},  and  thus  as  an  abstract  group  (sh  3s  z.  The  Fourier 


,_1.  i 2ninx 
(i>  ) = (e 


2 1 

decomposition  theorem  and  Plancharel  identity  takes  the  form  for  f e L (S  ) , 

1 

f = J i(n)f  (x)  , f(n)  = (f  ,f  ) = f ff  (x)dx  , 
ne'z  n n e " 


1. 


i.e.  f(‘)  is  a function  on  Z = (S  ) and 


(2.16) 


(f,f)  = / |f|2dx  = l |f (n) |2  = (f,f)  , 

0 nf  Z 


which  moreover  implies  f e L (Z) . 

In  addition,  since  (f  ) = np  , we  have  if  L = £ u . ) = P(tr^)»  that  Uf 

2tt  n n ] 2t  2n  r 

= P(n)  ' f , L(f)  = I f(n)P(n)^  (x)  , and  thus  by  the  previous  discussion  we  shall  think 
n n 

of  the  shift  operator  £ , and  ) as  analogous  quantities  in  the  two  Fourier 

2tt 

theories . 

The  discussion  of  the  Fourier  theory  of  Z has  already  taken  place,  for  from  the 
above  we  have  (S1)  = Z,  and  conversely  Z = (S1)  = S1,  for  as  before,  we 
may  interpret  the  characters  ^(-)  °f  S'*-,  as  characters  f ^ (x)  of  Z,  i.e. 
functions  of  Z.  Thus  if  we  read  the  Fourier  series  decomposition  and  Plancharel 

identity  for  backwards  in  (2.16),  we  have  the  corresponding  decomposition  and 

2 

identity  theorems  for  L (Z) . The  same  remarks  concerning  the  analogy  between 


-iD 

2n 


£,  ~XT  apply  < for  if  P (£ ) = \ a3  • £ -1 , then 

j£Z 

1 

P(£)^  = P(a  )f  , a = e T'1X,  P(£)f(n)  = / f p(a  )f  (x)dx  , 

n x n x J x n 

which  is  in  analogy  to  the  corresponding  formulas  (2.15)  for  both  Z involving  £ , 

n 

and  S1  involving  — . We  remark  that  its  sugqestive  to  "identity"  — — and 
2n  2 

j (£  - Z ),  but  that  identification  isn't  the  correct  one  for  our  purposes.  We  also 
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remark  that  the  only  definition  of  trace  for  an  operator  P (F  ) = £ ■ £ , is 

j6Z 

V 04 

tr  P = > (a  ) . , depending  on  which  case  we  are  in.  This  is  of  course  consistent 


3 


with  the  matrix  terminology  used  in  this  section.  An  important  question  answered  in 
the  next  section  is  the  corresponding  quantity  for  formal  pseudo-differential  operators. 
Remark  2.  We  note  that  formula  (2.7)  in  fact  defines  a natural  Poisson  bracket  for 

* 

•£  , as  is  well  known.  This  is  easily  seen  by  using  the  symmetry  of  < , >,  the  symmetry 


of  the  Hessian  H , an<^  the  Jacobi  identify  for  the  bracket  [ , ] . 


However,  given  a Poisson  bracket,  { , },  i.e.  a skew  symmetric  bilinear  derivation 
satisfying  the  Jacobi  identity  on  a C manifold  M,  one  automatically  has  a ‘stratifica- 
tion’ of  M by  symplectic  leaves,  i.e.  a ‘symplectic  stratification’.  By  stratification 
of  M,  we  mean  that  at  every  point  p e M,  there  exists  a smooth  leaf  £ 


(possibly  zero  dimensional)  running  through  p,  £ , and  moreover  £ is  part  of  a 

P P 


foilation  near  p which  foilates  nearby  £s’s,  perhaps  trivially,  and  in 
addition  M is  the  direct  union  of  these  leaves.  By 

symplectic  stratification,  we  mean  that  the  leafs  are  assumed  to  be  C°°  symplectic 


manifolds,  (or  points),  with  the  symplectic  two-form  u = to(p)  on  £ also  varying 

P 


smoothly  across  the  leafs.  In  other  words,  for  all  p,  given  any  smooth  vector  fields 
X,  Y defined  on  a neighborhood  of  p on  N , such  that  X(s),Y(s)  C T £ (s) , s e N 


P s 

(such  vector  fields  necessarily  exist),  u (X (s) , Y (s) ) (s)  is  smooth  in  s. 


To  show  that  we  have  the  above  'symplectic  stratification',  first  observe  that 
since  { f , g } is  a derivation  in  its  arguments,  it  is  local,  i.e.  its  value  of  p just 
depends  on  the  first  derivatives  of  its  arguments  at  p,  as  is  seen  by  substituting 
into  {.,f}  a function  which  vanishes  quadratically  fast  at  the  point  p,  and  using 
the  product  rule  of  derivations.  Next  define  the  Hamiltonian  vector  field  X,(g)  = {g,f}. 


and  since  we  may  restate  the  Jacobi  identity  for  { , } as  [X  ,X  ] = -X,-  - i < we 

ri  r2  ’ rl ’ r2 ’ 


see  that  A = UA  , A = span(X  ) 


P P 
peM 


f (p) 


, is  in  fact  an  integrable  distribution 


fee  (M) 

in  the  sense  of  Frobenius.  By  this  we  mean  that  near  p,  one  has  an  integrable 


J 
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distribution  fl  , s e N , such  that  ft  = A , I)  C t and  dimension  ( - ) = diraen- 
sp  ppss  s 

sion  (A  ) = n . This  is  clearly  the  case,  for  if  X^  (p)  ,X^  (p),...,X  (p)  span 

1 2 "P 

A , just  define  ft  = span(X  , ),  s « N . By  the  formula 

P s r , \ S ; TVS;  p 

1 n 

-1  fcl  t2  * t *i 

d<(iX  <Ji  = x , where  <J>  = <j>  »<(>  ....  • <J>  n,  <{>  the  flow  generated  by  a vector 

fo* 

field  of  the  form  X^,  we  see  that  a leaf  cannot  change  its  dimensionality,  which  proves 

the  assertion  that  M equals  the  direct  sum  of  the  leaves.  We  call  the  leaf  running 

through  p,  E . We  define  the  two  form  on  E^  at  p by  w(Z,Y)  = {f,g}(p),  where 

Z = Xf  (p)  , Y = X^(p).  Clearly  u is  well  defined.  To  see  dco  = 0,  just  apply 

(using  the  intrinsic  definition  of  d)  , X (p)  , X (p)  , X (p)  to  dco,  and  use  the 

2 1 2 3 

Jacobi  identity  coupled  with  d f = 0.  That  w is  nondegenerate  is  automatic,  since 

Tl^(p)  = A . This  proves  our  assertion.  Note  that  given  a 'symplectic  stratification' 

of  M,  one  automatically  has  a Poisson  bracket,  for  simply  define 

{f,g}(p)  = {f,  ,g.  },  , where  f,  = f restricted  to  A , etc.,  and  { , } 

1 A !A  ' A p An 

p p p p * 

denotes  the  Poisson  bracket  A^  inherits  as  a symplectic  manifold.  That  { , } 

is  a Poisson  bracket  follows  from  { , being  one,  and  the  stipulation  placed  upon 

P 

the  smoothness  with  which  the  symplectic  two  form  on  the  leaves  may  vary. 

We  also  note  that  if  we  assume,  for  instance,  that  M is  not  compact  and  X^  = 0 
0° 

implies  M < Cq(M)  , the  functions  of  compact  support,  then  if  G is  the  group  with  Lie 

. oo  * 

algebra  A = {X^.|f  e Cg(M)},  we  may  think  of  G x M C t G,  via  the  pairing  <p,X^>  = f ( p)  , 
and  the  above  construction  as  an  example  of  the  Kostant-Kir illov  orbit  construction. 

This  will  be  discussed  elsewhere. 

Thus  given  a Poisson  bracket,  it  is  natural  to  compute  the  symplectic  leaves  it 
induces.  In  the  case  of  (2.7),  the  leaves  are  nothing  more  than  the  orbits  of  Kostant- 
Kirillov.  It  will,  however,  not  always  be  the  case  that  a finite  dimensional  group 


generates  the  orbit  structure  on  M. 


3 . The  Generalized  Kor fcowcg-deVr ies  Equations 

In  order  to  apply  the  considerations  discussed  in  Section  2 to  the  generalized 
Kor teweg -'deVries  equation,  we  need  to  define  the  appropriate  formal  Lie  group  G.  To 


this  end  we  introduce  the  commutative  ring  R = R(..,a 


,art, . . ,a  . , . . ) = R(a)  over 

o 3 


the  complex  numbers  equipped  with  a natural  derivation  D,  consisting  of  polynomials 
in  the  a^  and  their  derivatives.  One  defines  the  'indefinite  integrals',  I = R/DR, 
i.e.  R modulo  DR,  and  we  shall  use  = for  the  equality  sign  in  T.  Note  that  D 
treats  I elements  as  if  they  were  constants. 

We  now  define  the  ring  of  formal  asymptotic  pseudo-differential  operators  over  R, $, 
to  be  the  ring  of  formal  Laurant  series  in  the  variable  £,  where  we  may  heur istical ly 
think  of  £ as  contained  in  a real  neighbourhood  of  <»  which  is  moreover  embedded  in 
a complex  neighbourhood  of  °°,  and  we  shall  take  the  coefficients  of  the  series  to 
be  R elements,  and  hence  D extends  to  4>  naturally.  We  shall  also  require  that 
each  series  terminates  at  some  £ ^ , j < <».  The  ring  multiplication  shall  be  just  th»  ont 
used  in  pseudo-differential  operator  theory,  namely 


(3.1) 

Note  we  may  write 


= l • (-iD)V<t2  • 


v>0 


(3.2) 


<f  = { 4>  = J a^£1|a,  a generator  of  R(a)}. 


— <»<  i<N<°o 

The  motivation  for  these  definitions  was  discussed  in  the  introduction,  and  we 


reiterate  that  we  heuristically  think  of 


elements  as  germs  at  °°.  We  briefly  discuss 
d 


(3.1).  If  R elements  were  just  C functions  of  x,  D = - — , then  the  differential 

dx 

ixf  i x r ixT 

operator  P(-iD),  with  P a polynomial,  acts  on  e ^ via  P(-iD)(e  *)  = P(£)  • e 


and  clearly  (P^^  (-iD)  (P2  (-iD)  (e  )))  = (P^O  » P2  (£  ) ) 


ix£ 


with  o defined  in 

(3.1).  Thus  $ models  the  image  of  the  algebraic  isomorphism  implicit  in  the  above 
comments,  P(-iD)  H-  P(£),  in  a formal  way,  and  moreover  extends  the  image  of  the 
isomorphism.  These  comments  also  explain  why  o is  associative. 


t 


The  care  taken  in  the  definition  of  ♦ was  in  order  to  define  the  trace  functional 
< > : 4 ►*  I , by,  assuming  <p  = £ a.^1, 

(3.3}  tr  ^ =((})=  a ^ e I, 

where  if  <j>  e R,  <(>  t I denotes  the  element  in  I with  representative  element  ip  . We 
now  state  the  simple  but  fundamental  result  concerning  < >. 

Theorem  2.  If  = ^ ° ♦ 2 ” ^ ° ^ 1 ' we  ^ave 

(3.4)  ( > = 0 . 


Proof : It’s  here  that  the  duality  between  x,E,  in  Fourier  theory  plays  a 

crucial  role,  as  evidenced  by  its  appearance  in  the  role  of  multiplication  (3.1).  Its 
immediate  from  (3.1)  that 


4>2 1 


hence  pro jecting into  I = R/DR  via 


a + DB,  a, 8 e 4>  , 


3C 

we  have 


[V*2]  = If  ' 
.-1 


and  so  clearly  can  have  no  £ term,  which,  by  (3.3),  concludes  the  proof 

of  Theorem  2. 

We  also  remark  that  one  can  prove  Theorem  2 by  direct  computation,  which  is  useful 
for  Section  4 but  then  one  fails  to  see  why  our  choice  for  ( > is  really  unique. 

As  a consequence  of  (3.4)  we  make  the  important  definition,  analogous  to  (2.1)  of 
( , ) , 


(3.5) 


(A  , A ) z(d) 

vl’v2  V1 


,>  - 


and  so  by  Theorem  2,  we  have  that  ( , > is  symmetric  in  its  arguments,  which  is  the 

crucial  point  in  order  that  the  operation  [<(>,•]  be  skew-symmetric  with  respect  to 
( , >,  which  played  an  important  role  in  the  considerations  of  Section  2. 

We  now  define  ♦ to  be  the  ring  over  4 of  formal  Laurant  series  in  which 

we  shall  heuristically  think  of  as  being  a member  of  the  complex  numbers  in  a deleted 
neighbourhood  of  °°,  minus  the  real  numbers.  4>^  inherits  its  multiplication  structure 
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from  *>  and  the  law  of  exponents.  We  now  single  out  # , a submanifold  of  <J>, 


defined  bv 


$ = { <t>  - r + 7 a . £ a a generator  of  R(-)l  . 

n . „ 1 1 i 

0*1 <n- 2 


We  can  now  define  the  resolvent  operator  6? 


by  requiring 


<W  : $ h-4»,  x = A - £ , (a)  a,  , 

n n r n *A 


(A  ~ <£)  - (A  - 4> ) o $ ^ - 1 


We  compute  in  the  following  way,  writing  $ « £ * 


(A  - <t>) 


{ (A  - £ ) - <t>]  = 4>  • (A-rn)  - 0^  0 4>=1 


and  so  satisfies  an  equation  of  the  form  (I  - T)x  = y,  with  x = <fe  ^ , y = (A  - £°)  * , 


(A  - £ ) 


, hence  $ = £ T^(A  - £n)  \ from  which  it's  easy  to  see  by  (3.1)  that 


(3.8)  1 = £ * . (-1) " Cm(?n  - A) 

A m,  Jc 


_ _ m + £ 

m > 0,  i > 0,  an  integer  , 


for  <t  . € R - R (a  .a  ...  a ),  where  R is  the  ring  whose  elements  are  finite 
m,  x,  n n 0 1 n-2  n 

sums  and  products  of  a^,  0 i <_  n - 2,  (see  (3.6)),  and  their  derivatives.  In  (3.8) 

we  use  the  notation  of  Gelf and-Dikii  [1].  Using  the  construction  of  (3.8),  we  can  mimic 

the  construction  of  R.  Seeley  [7],  in  a formal  way,  namely  we  define  for  'nice'  f, 

(3.9)  f(<f>)  = [f  (A  )<f>  ] 

* (A,5n) 


where  for  h = h(X),  [h] 


is  the  formal  residue  term  of  h at  £ , i.e.  the 


coefficient  of  the  (A  - £ ) term  of  the  Laurent  series  in  (X  - £ ) , given  by 
formally  expanding  h = h ( A ) about  £n.  We  keep  this  notation.  This  of  course 

necessitates  expanding  f(x)  about  x = £n,  which  places  restrictions  on  the  choice 
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of  f.  One  uses  (3.8)  to  compute  (3.9).  In  practice  one  may  have  to  extend 

l1  I1  to  include  flip).  As  a simple  example,  by  (3.8),  (3-9),  one  easily  computes 


E * 


where  A 1 P , I 1 - For  the  case  that  s = - , N e Z,  but  not  a multiple 

m 

of  n,  we  have 


N 


N 


(3.10) 


f " = 0 " £>Q  Vn^'  h6nCe  = Vl,N  ‘ ° 


where  we  may  heuristically  interpret  o 


f sign(C),  n-odd 
VI  , n-evenj 


but  it  is  really  just 


an  adjoined  symbol  obeying  the  formal  rule  of  multiplication  a =1. 

N 

We  amplify  (3.10)  . Here  we  think  of  cjn  € and  extension  of  $»,  $ = $ e a * <j. 

We  must  think  of  a as  just  an  adjoined  symbol,  obeying  = 1,  but  if  we 

interpret  ) f $ as  a series  in  real  £ embedded  in  a neighbourhood  of  °°  in  the 

complex  plane,  then  o has  the  interpretation  given  above.  Implicitly  we  have  modified 

the  ring  R,  i.e.  R -*  R * o • R,  hence  I -<•  I ® o • I.  That  the  trace  has  all  of 

its  usual  properties  follows  from  ^-o  = 0.  In  general  <f> S e £ e , with 

T,sf  C 

°x  * °_  = °T  . , (C  thecomplex  numbers).  Hence  it’s  obvious  that  if  s is  irrational,  *S 

1 2 1 2 

has  no  residue  term  and  so  < * S > = 0,  s * N f Z.  In  the  future  it  shall  be  understood 

we  are  working  with  the  enlarged  ring,  R ® a • R,  but  we  shall  not  change  notation. 

In  addition,  we  define  K(4>,X),  using  the  notation  of  (3.9),  and  (3.8),  as 


(3.11) 


K(<M)  5 [♦,]  1/n  i l u.<xV(n+ll"1)  fi  . u e I . 

* <£.A1/n)  l>0  * (X1/n)  * 


This  is  easily  seen  by  computing  the  formal  residue  of  terms  of  the  form 

m n n 

Tm,n  = 5 in  (3.8),  which  equivalently  comes  down  to  expanding  t (£) 

' m,  n 


-16- 


in  powers  of  ( £ - A*^n)  and  picking  off  the  (£  - A*^n)  * coefficient;  in  other 


words,  we  think  in  terms  of  expanding  = 4^(S)  about  A 


1/n 


Notice  we  have  formally 


declared  (A^njn  = A in  t rather  than  "pick  a branch".  In  analogy  to  Section  2, 

m,  n 

we  define  A.  . = { Y cr},  and  the  natural  projections  P,  . into  the  A.  . 's. 

i»3  sZ:i  s,Ys  J x,3  i,j 

We  shall  work  with  the  formal  Lie  group  G = 1 + A ^ , with  formal  Lie  algebra 

£ = A __  Alternately  we  shall  choose  to  represent  jC  in  the  following  form 

(3.12)  £ = { l (e  - iD)"j'V  1$.  e R}  , 


j>0 


3 3 


where  (£  - iD)  3 \|> 


r -i-l-v  i+V  V 

l C (iD)  <\i,  j ^ 0.  It  is  obvious  that  the  two  defini- 


v>0 


tions  of  £,  are  the  same,  in  fact  if 
(3.13) 


I - l (£  - iw-J-V  - l l^ldDlVe'1 

>0  j>0  J v>0  J 

j>0 

= i < ! , 


-l-3-v 


k.  k 

(3.14)  c.  = 7 k (iD)Vb,  , and  so  b,  = 7 |k)(-iD)Vc 

k v k-v  k . . I v / 


k-v 


k>0  v=0 

we  must  have 

k 

I ( iD)  Vb,  , and  so  b,  = 

_ 1 v / k-v  k , . 

v=0  k=0 

The  second  relation  in  (3.14)  being  an  easy  consequence  of  the  first.  We  shall  denote 

j 

I 

s=k 

A . . We  identify  the  dual  of  £ , 

k,3 


by  A - { £ (€  - iD)  • <(>  1 4>  e R} , and  by  P . the  natural  projections  into 

k , 3 e=v  s s k , 3 

£, 

(3.15)  £ = { l d> . e1 1*  . e R}  = A 


oo>N>i>Q 


If  A = l a 

i>0 


iC1  e i , B = l (5  - iD) d . 

j>0  -1 


0," 


. , then 


(3.16) 


< A,B>  = <A°B)  = 7 a , b . 

• „ 1 l 

i>0 


as  easily  follows  from  (3.15),  (3.1).  Note  this  is  in  complete  analogy  with  the  formula 
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following  (2.10).  In  other  words,  (E  - iD)-]  1 is  due  to  £-,  j ^ 0 under  < , >,  i.e. 

(3.17)  (£?>*  = ( (£  - iD)  -3’1],  ( (E  - iD)  = £ ? , j > 0 ; 

where  * denotes  duality  under  < , >.  Heuristical ly  'total  integration'  j + 3 tiroes 

l 

is  dual  to  'differentiation'  j tiroes.  We  now  are  in  a position  to  state  the  main 
theorem  of  this  paper,  but  we  first  mention  that  we  shall  think  of  covariant  tensors  as 
multilinear  operators  on  vector  fields,  linear  over  I, 

Theorem  3.  If  A e 0^  = ( [g  ^Bg)+|g  e G},  where  + denotes  projection  onto  Aq  < 

and  B e An  ^ i.e.  00  is  just  the  orbit  of  B under  the  dual  of  the  adjoint 

m 

action,  and  so  A = £ a.£  , then  the  (formal)  natural  symplectic  2-form  w on  0 

i-0  1 B 

is  given  by,  at  the  point  A,  with  l^e  £, 

w([A#£1]+, [A,*2]+)  = < A, [4  ,£  ] > = < [A,*  ]+,£  > . 

Thus  Hamilton's  equations  associated  with  <*>,  for  the  Hamiltonian  H = H(A)  e I,  a 
polynomial  in  a^,  i = 0,1,..., m and  its  derivatives,  is  given  by 

(3.18)  A = XH  = [HA,A]m_2  , 

( ]n  denoting  the  natural  projection  onto  , where 

0,  n 

m 

(3.18)  Ha  = I (5  - iD)'3'1  , 

3=0  3 

DH  „ 

with  being  the  formal  variational  derivative  of  H € L with  respect  to  A..  In 

3 3 

addition,  the  (formal)  Poisson  bracket  { , },  is  given  by 

(3.20)  {H,F}  i < A, [H  ,F  ] > . 

Finally  this  Poisson  structure  agrees  with  the  structure  of  Gelfand-Dikii  [1]  if  we 

take  a = 1,  a , = 0.  We  can  do  this  as  a , a , are  orbit  invariants.  At  this 
m m-1  m m-1 

point  we  may  restrict  our  differential  ring  R to  only  consist  of  the  polynomials  in 

the  a^'s,  i = 0,1, ...,ra,  and  their  derivatives,  over  the  complex  numbers,  with  the 

adjointed  branch  symbol  a,  i.e.  R = R(a  ,a, ,...,a  ). 

0 1 m 

Proof : We  are  formally  in  the  same  situation  as  in  Section  2,  and  so  we  shall 
not  repeat  the  discussion  now  but  merely  amplify  certain  points,  and  hence  just  sketch 
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the  proof.  For  instance  if  H = H (A)  r I,  tiie  gradient  of  H with  respect  to 

( # > is  first  defined  uniquely  by  requiring 

X (H)  = < H,X>,  H e JC  , 

A A 

r D * 

where  X = l — — , is  a vector  field  on  £ , which  in  the  inner  product  is 

i>0  i 

* r'  i 

naturally  identified  with  the  £ element  \ X . f,  . For  H depending  only  on 


a0'  ai ' * " * 'am'  We  automatically  have  e A j#  and  then  formula  (3.19)  follows 

from  (3.16).  Also  see  (4]  for  amplification  on  related  matters.  The  formulas  for  6 ,io# 

B 

Hamilton's  equations,  { , },  were  clearly  taken  over  directly  from  the  discussion 
of  Section  2. 

We  emphasize  that  these  equations  can  only  be  taken  in  a formal  sense,  which  implies 
that  u>  is  formally  closed,  (using  the  intrinsic  definition  of  dco)  and  therefore  that 
{ , ,l  is  a Poisson  bracket,  i.e.  is  a skew  symmetric  derivation  which  satisfies  the  Jacobi 
identity.  We  amplify  this  point,  taking  the  point  of  view  of  the  formal  calculus  of  variation  of 
Gelfand-Dikii  [see  4] . We  define  the  operation  of  exterior  differentiation,  d,  using 
the  intrinsic  definition  which  relies  on  the  Lie  bracket.  This  definition  gives  us 
no  trouble,  as  our  ring  R(a)  is  composed  of  polynomials.  For  instance  for  H € I, 
dH  = E da..  From  the  formula  for  go  given  above,  one  verifies  directly  that  dco  = 0, 

as  indeed  is  true  in  the  case  of  a bonified  Lie  group.  In  the  usual  verification  that 

{H,F}  = to(XH,Xp)  satisfies  the  Jacobi  identity,  one  used  that  doj  = 0,  and  that 
2 2 

du>(HH,-)  = d(dH)  = d H = 0.  That  d =0  follows  from  the  intrinsic  definition  of  d, 
and  says  in  this  case,  in  coordinates,  that  the  2-form 


d(dH)C,-)  =d(S^-da.,  = l To  (•)  - (•) 

1 i<3  L 1 3 J 1 J 


A da,  =0  , 


t j 

where  <p , <P  c R,  is  the  (formal)  directional  derivative  of  <J>  with  respect  to  a.,  or 

i 

equivalently  d2H(Z,Y)  = < Z, 6H„  (Y) > - < Y, 6h  (z)  > = 0) , 6 = ~ , hence 

A A 0A 


& DH  6 DH  , . 

7—  = TT-  ^ — (*)  , all  l,]  . 

oa  . Da  . oa  . Da  . 
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This  is  just  the  well  known  fact  that  the  directional  derivative  of  is  a symmetric 

operator  with  respect  to  < , >.  See  [4,3]  for  amplification  on  this  particular  case 

of  the  Poincare  Lemma,  where  this  fact  is  discussed  and  proved  differently.  One  can 

also  show  directly  from  (3.20) , that  { , } is  a Poisson  bracket.  One  uses  for  this 

the  symmetry  of  < , >,  the  Jacobi  identity  for  [ , 1,  and  that  the  formal  Gateaux 

derivative  of  is  a symmetric  operator  with  respect  to  ( , >,  which  is  precisely 

what  one  would  use  in  Section  2,  as  was  mentioned  in  Remark  2. 

Finally  the  remark  concerning  the  Poisson  brackets  of  Gelfand-Dikii  follows  first 

from  the  fact  that  since  a , a , are  orbit  invariants  they  may  be  taken  equal  to 

m m- 1 

1#  0 respectively.  Note  that  in  general  the  symplectic  structure  will  depend  on  a , 

m 

as  indeed  the  following  computation  will  make  clear.  We  now  compute  (3.18) 
using  formula  (3.1),  which  will  finish  the  proof  of  Theorem  3.  First  observe,  setting 

-f-  = H 
Da . a . 

3 3 

»„  - I ,t  - - I , 

j>0  3 3>0  ] 


A = IHa.A] 


-j-l-sf  j+sl  s 


(iD)  Ha  ) « (1  ak£  ) 


a akek)  - a e'j_1-sp 


-s  j:si(iD)sH 


(iD)SH  (iD)Va, 
3 k 


. - j-l-s-v+k 


(with  -j-l-s-v  + k = r^0,  v+s  = y) 


l :l' ■ «*»"-< 

3 


r+ j+u+1 


_ r+j+u+ll  j+sl,  v , . u 

E v s ('X)  ar+j+u+l(lD) 


-i+W)  (iD) W (H  • a ) 

U / a^  r+j+u+1 


r-0  /r+u 


+lJ  a , . (-iD)  ^H. 

V I r+3+p+l  a^ 


Cr=  I 


a/. 


/ 


where  the  sum  in  the  brackets  is  understood  to  extend  over  ull  nonnegative  integers  with 

the  proviso  that  only  terms  with  a_.,  j = 0,1,..., m,  can  appear.  This  formula 

agrees  with  the  formula  of  Gelf and-Dikii  11]  if  we  set  a =l,a  , - 0,  i.e.  in  the 

m m-1 

setting  of  (3.6),  and  so  Theorem  3 is  proven. 

Remark  3.  In  fact  in  [1],  Gelfand-Dikii  arrived  at  this  formula  from  computational 
considerations  which  shall  become  apparent  in  Theorem  4,  and  then  proceeded  to  verify 
by  a difficult  computation  that  (3.21)  in  fact  yielded  a Poisson  bracket  via 

(g,h)  = ^7 1 , with  H the  Hamiltonian  in  (3.21)  . 
at,H 

Note  (3.21)  is  equivalent  to  stating  that  the  Hamiltonian  vector  field  associated 


with  the  function  H is 


, D [ D D 

where  — = - — , , with 

Da  Da  Da 

1 0 m-2 


J = A - A , 


where  A is  the  {m  - 1)  * (m  - 1)  matrix  differential  operator  with  components 


A -If 

rs  u 


i ay+l+r+s('iD) Y'  0 < r < m - 2,  0 < s < m - 2 , 


and  A is  the  dual  of  A with  respect  to  the  natural  inner  product  on  R m , 

n 

m-2  

(v,w)  i v.w.  6 I,  and  the  remark  after  (3.21)  concerning  summation  applies  above. 

Clearly  for  m = 2,  X = -2iD  • , the  well-known  Poisson  structure  of  Gardner 

H Da0 

associated  with  the  Korteweg-deVries  equation.  Equivalently  the  Poisson  bracket 


{g,h}  = xhg=  (A5H  _(5H 

H 1 Da  Da  l Da 


Remark  4 . It  is  interesting  to  study  the  orbit  of  £m  + £ a.^1  under  the 

i=0  1 

group  G,  in  fact  we  claim  that  a , a are  always  orbit  invariants.  Since  a 

m-2  m-3 

typical  element  of  G is  of  the  form  1 + x,  x f A ,,  i.e.  x = T x f”3.  to 

-"(-l  . , i 

g>i 

prove  the  remark  we  compute,  (using  (3.1)),  that  under  the  action  of  G 


U + X + X2  + •••}) 


[ 


6A  = ((1  - x) A ( 1 + x + x + •••)!  - A = [[A,x]+  - 

“ (c1Dx1)f,m_2  ♦ (c2x1dxi  + c3D2Xl  + c4Dx2)Cm-3  + n.  n e A0tn^ 


with  C,  , C . C,,  C.,  constants  depending  on  ra,  and  so  6a  _,  6a  , C DR, 
1 2 3 4 r m-2  ra-3 


s 

thus  proving  the  assertion.  Tn  fact  more  is  true,  namely  (A™),  1 £ s < m - 1 are 
also  orbit  invariants,  which  in  fact  implies  the  above  assertion.  (See  Corollary  1 at 

the  end  of  this  section.)  Thus  in  the  case  of  the  Korteweg-deVries  equation,  an  orbit 
2 

of  -D  + q is  specified  by  q = constant.  Thus  the  apparent  one-dimensional  degeneracy 
of  the  Gardner  structure  is  really  a consequence  of  the  representation  of  an  orbit  by 
q,  which  really  contains  a hidden  relation,  q = constant.  This  phenomena  is  well 
known  in  the  Toda  system,  where  as  previously  discussed  in  Section  2,  one  has  the 
similar  orbit  relation,  trace  A = constant. 

Remark  5.  Suppose  we  had  instead  defined 


* 

n 


m-2 

+ £ (a.f,3  + (f  - iD)3a.)|a.  a generator  of  R(a)),  m < 

j=0  3 33 


i.e.  we  coordinize  ^ by  these  a^'s.  Then  as  our  arguments  have  shown  (3.18) 

defines  a Hamiltonian  structure.  Of  course  we  have  to  compute  H , as  in  (3.19),  or 

equivalently,  we  must  compute  the  dual  of  the  = 5 3 + (£  - iD)3,  0 < j < °°  in  .£  . 

The  computation  of  the  dual  of  £3  depended  on  inverting  the  first  relation  in  (3.14), 

but  in  fact,  it  was  only  necessary  to  invert  the  linear,  strictly  triangular  relation 

up  to  degree  m + 1 to  compute  the  Hamiltonian  structure  from  (3.18).  This  is  always 

a trivial  matter,  since  m is  finite,  and  the  relation  to  be  inverted  is  upper 

triangular.  Hence  we  can  use  the  above  remarks  to  compute  the  Hamiltonian  structure 

of  in  the  coordinates  of  the  'formally  self-adjoint  operators'  i.e.  with  the  above 

coordinizatior,  of  ♦ . In  fact,  all  the  results  of  this  paper  have  and  shall  be 
m 

independent  of  which  choice  of  coordinates  one  uses  on  $ , although  for  clarity  we 

m 

shall  sometimes  just  work  with  the  original  coordinates  on  <J>  . However,  in  the  case 

m 

of  the  formally  self-adjoint  operators,  we  have  extra  reality  conditions,  if  we  require 
the  aj's  to  be  real,  for  instance  the  Poisson  bracket  and  traces,  (3.10)  must  be 
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real.  Hence  for  the  rest  of  the  paper,  we  shall  mean  either  of  th*  two  coord inati zat ions 
by  out  expression  unless  otherwise  stated.  For  the  * sol f-ad joint  * rase,  we  have 


for  J,  the  symptectic  matrix  operator,  for  m = 3, 


0,P 

J = 

D , 0 ‘ 


which  appl  .-s  to  the 


Boussinesq  equation,  see  [5]  , [f»]  . For  m = 4 , one  can  show 


J = 0 -1/2D  0 , k = -1/2D  f a.D  f Da0,  see  [5]  . 

0 0 0 

L. 

Remark_fi.  If  we  had  instead  chosen  R * R = R (.  . ,3^  j)  , . . ) = ML(n,R>  , the  ring  of  n*n 

matrices  over  R,  instead  of  R to  tie  the  base  ring  of  our  construction,  then  we  would 

have  to  redefine  for  4 e 4,  (see  (3.3)),  tr  $ = matrix  trace  a t I,  which  extends 

definition  (3.3).  Note  r is  still  R/DR.  All  the  constructions  and  arguments  of 

k * 

this  paper  apply  to  this  case.  For  instance  (E.  .£ . ) , where  E . is  the  matrix  with 

13  13 
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its  1,3  component  unity,  all  other  elements  zero,  equals  ('  - i D)  . E , and  1 

ji 

shall  often  denote  the  n * n-identity  matrix,  hence  H is  the  matrix  such  that 

ak 

DH 

[H  ] . 0 = - — q-  which  would  have  to  be  substituted  in  (3.21).  In  addition,  one 

ak  11 1 Da/0) 

k 

may  specialize  the  coefficients  of  A to  lie  in  some  ring  of  matrices,  which  amounts 
to  changing  the  base  ring  to  that  ring  of  matrices. 

We  now  pursue  the  result  analogous  to  Theorem  1,  but  first  some  brief  observations. 
Suppose  H = <AV>,  v = ^ , N an  integer.  We  wish  to  compute  H ^ . We  first  write 


A = [X  *X1 


, by  (3 . 9) , and  so  ) . In  addition,  by  (3.7)  we 

a,cn>  x lu,cn) 


have  = 4 o 6A 


X'  dA  “ n 


= “V5,  , with  6 indicating  an  increment,  and  so 


<5h  i < (XV4.  ° <SA  . 4]  > i ( [xV]  , , 6A  > + < 3_  { • } ) 

X X (X,5n)  X l(X,Cn)  5 

- < [aV]  , a>  = < [-  (4r)xV]  -SA> 
x (X,f ")  L 3>  ' J(X,£n) 

[v  _ 

trr~  •4,1  + 1 3 . ( ' 1 ] , 6 A > = < [ vA  V ^4  ) ,6A  >=<  vAV  1 , 6A  > 

3X  XJ(X,^)  A (X,5n)  x (X,?n) 


(when  we  have  repeatedly  used  that  a perfect  derivative  can  have  no  residue  term) , hence 


we  have 


<AV>  = (vAV"1)  , 

A - (m- 1 ) 


A e $ , 


where  the  operation  ( ) denotes  the  natural  projection  into  A , by  P 

~n  -n,  -1  7 -n , 0 

We  next  compute  with  v = H ==  < AV } = < [xV]  ) i (interchanging  residue 

* (X,£m) 

evaluations,  see  (3.3),  ~ [ X VM<t> , A , ] = a ■ by 

(c  / A ) 

(3.11).  In  this  calculation  we  used  the  meaningful  convention  that  the  residue  term 
at  00  is  -A  . One  can  justify  the  interchange  of  residue  computations  in  the  above 
by  thinking  of  performing  the  above  in  some  analytic  setting  under  appropriate  restric- 
tions, where  the  residue  calculations  are  just  achieved  by  contour  integrations  in 
A,£  in  a small  sector  near  00 , and  the  results  in  this  setting  imply  formal  identities 
which  yields  the  above  formal  result.  We  have  thus  shown 


K W>,X)  = a ■ l ( A m > (Xm  ) 


e-i  1 

m . m - (m+£-l,  . 


a • < an*,  - (* ) n_1>. 


<3.24,  <A  m > - l^](Am,  , n . 

A l m / - (m-1) 

Note  that  certainly  the  computational  equivalents  of  Lemma  1 can  be  found  in 
Gelfand-Dikii , [1]  but  of  course  proved  and  viewed  in  a different  fashion.  Note  that 
(3.23)  implies  that  the  A , j = 1,2,...  terms  of  K($,A)  are  missing,  which  is  no 

surprise.  We  can  now  prove  the  analog  of  Theorem  1,  namely: 

N+m 

Theorem  4.  (Gelfand-Dikii) . If  H = H = — — — — ( A m ) , A e $ then  Hamiltons 

N + m m 

equations,  A = [HA,A]m  ^ , imply  the  Lax-isospectral  equation 

N 

(3.25,  A = [A.P  I,  with  P = P„  = (a”, 

N N + 

with  ( , as  usual  indicating  projection  onto  . The  H is  — 2 — times  the 

+ 0,«°  N N + n 

1^ 

coefficient  of  the  (Xm  , term  of  the  'trace'  of  the  resolvent,  K(<f,X).  In 
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addition  the  ( A ) are  in  involution  with  respect  to  ( , },  the  Poisson  bracket. 

Proof : Theorem  4,  like  Theorem  1,  depends  on  the  crucial  observation  (AS,A1  - 0 

for  all  s.  This  is  a consequence  of  the  definition  of  f(L)  by  residues  (3.9)  and 

the  well-known  functional  equation  of  the  resolvent  — = <p  o <j  , which 

A1  " A 2 A1  A2 

N N N 

implies  A * A = A^  . Clearly  A a (A  ) + (Am)  , with  ( ) denoting  projection 

N N N N N 

into  A-oo,-i'  and  so  0 = [Am,A]  = [(Am)+  + (Am)_,A],  therefore  I(Am)_,A]  = [A,  (a"1)  +] 
and  so  we  have 


(3.26)  [A,  (A  ) ] - [ (A  ) . /A)  + [0  ,A]  , 

+ - (m-1 ) m 

where  0 is  contained  in  A .We  project  (3.26)  onto  A n , observing  that 
-m  -°° , -m  0 , 00 

N N 

[A#  (A  )+]  € while  e [(A  e A m_2*  We  t^lus  after 

projection  onto  A ^ , arrive  at 

0 , 00 

N N 

(A, (A  ) ] = [ (Am)  ,A]m  2,  which  is  a member  of  _ , 

+ -(m-1)  0,m-2 

and  thus  as  a consequence  of  Lemma  1,  Theorem  4 is  proven,  modulo  the  involution  statement. 

We  now  prove  the  involution  statement,  namely  that 

N K 

{(A  >,<Am>}  = 0,  N , K integers  , 

where  { , } is  the  Poisson  bracket,  (3.20),  i.e.  we  show  that  the  traces  of  A are  in 
K N K N 

involution.  If  H = (A  ),  then  {(A  ),(A  )}  = ( A ) , by  the  Hamiltonian  formalism, 

K 

with  H = < An>,  i.e. 


{<  A >,(A  >}  - <A,<A  ) , by  the  definition  of  gradient  , 


- < [A,PK1 , (Am)_(m_l) > = < [A,PK] , (Am)_>,  (since  [A,PK]  e AQ  m_2>  see  (3.16)) 


r 

t 


S 


N N 

- < [A,P^]  (A01)  > = < [ A, P^ ] Am > , (since  the  trace  only  picks  out  the  £ 

N 

= ( P , [A™, A]  ) = 0,  thus  completing  the  proof  of  Theorem  4. 


We  wish  to  point  out  the  parallels  between  the  proofs  of  Theorems  1,  4.  We  also 
wish  to  remind  the  reader  of  Remarks  5,  6 which  imply  that  Theorem  4 really  constitutes 
a generalization  of  the  theorem  of  Gel ' f and-Dikii , as  it  includes  the  cases  discussed 

in  these  remarks,  with  but  minor  modifications  in  terminology. 

s 

Corollary  1.  If  A e $ , < Am> , s = 1 , 2, . . . ,m  - 1 , are  orbit  invariants,  which 
m 

in  fact  implies  a _,  a . are  orbit  invariants. 
m-2  m-3 


N 

Proof:  Since  P„  = (A™)  , P.  = 0 for  N a negative  integer,  hence  [A,P  ] = 0, 

N + N N 

and  so  in  Theorem  4,  A = 0,  hence  A = [HN  ,A]m  2 = 0 for  N'  < m.  But 

A 

(Ha# A] m 2 = 0 is  a necessary  and  sufficient  condition  for  a function  H to  be  an 

orbit  invariant,  by  (3.20),  and  thus  we  have  shown  that  <An>,  s = l,2,...,m  - 1, 
are  orbit  invariants. 


m-2 

To  see  that  a a _ are  orbit  invariants,  observe  that  A = £ + Y a.E 

m-2  m-3  . l 

1=0 

has  weight  m if  we  assign  £ weight  1,  and  a.  weight  m - i.  From  the  definition 
N N 1 

A1",  (3.9),  it  is  easy  to  see  A™  has  weigh  N if  we  further  define  D^a.  to  have 

N 1 

weight  j + (m-i) , i.e.  D to  have  weight  one.  Thus  if  <Ara>  = E^,  must  have 

weight  N + 1,  and  so  we  must  have 

1 2 


<A  > =Ciam-2' 


< A‘">  = C a 

2 m-3 


C^,  being  (nonzero)  constants  which  can  be  computed  explicitly.  We  conjecture  that 

the  algebra  formed  by  the  invariants  of  the  corollary  provide  a complete  list  of  orbit 
invariants,  at  least  those  which  can  be  expressed  as  I elements. 


*3 


term) 
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4. 


Lenard  Relations 


We  are,  from  the  considerations  of  the  previous  two  sections,  in  a position  to 

compute  recursion  relations,  or  so-called  Lenard  relations  for  various  quantities  of 

interest,  sec*  [5]  for  instance,  where  such  relations  are  discussed. 

We  consider  the  setup  in  Section  3.  The  crucial  too]  in  these  relations  is  the 

identity  A A - A'  , which  plays  such  an  important  role  in  Theorems  1,  4 in  the 
s 1 

weaker  form  (A  ,A  J = 0,  and  whose  proof  was  discussed  in  the  proof  of  Theorem  4. 

In  addition  the  general  relation  < lE,P]  > = 0 is  important  (see  Theorem  2) . 

We  make  the  usual  natural  decomposit ion  of  $>  associated  with  £,£  , namely 


B € $ implies  B = B + B , with  B e £ = A , B e £ ' A 

+ “ + 0,00  - “00,-1 


• We  thus  have, 


with  A e $ 


however. 


(A  ) + = (A  • A)+  = (A  ) + * A + ((An)_  • A)  + 


(A  ) - (A  ) + 9 *=  H + 0 , 

-(n-1)  -n  A -n 

with  e , by  Lemma  1,  using  the  notation  of  Theorem  4.  The  above  yields, 

still  using  the  notation  of  Theorem  4,  (note  PN  = 0,  N a negative  integer) 

(4.1)  P = P A + ( * A ) + E , 

N+n  N A + 0 


with  Eq  equaling  the  coefficient  of  the  (£  - iD) _n  term  of  (An)  , i.e.  the 


coefficient  of  the  component  of  [ (A  ) - H 1 in  A .In  addition,  since 

A -n.-n 

N N 

( r(An)_,A]  > = 0 by  Theorem  2,  the  £ ^ coefficient  of  f (An)  ,A]  is  an  exact  deriva- 

N 

tive,  identically  in  the  coefficients  of  (An)  ,A.  But.  by  the  rule  of  multiplication, 
(3.1),  the  coefficient  of  the  £ ^ term  of  the  bracketed  expression  is  of  the  form  (up 
to  a constant  multiple) , 

N 

DE  + Q(( An)  ,.)  = DE.  + Q(h”) , by  Lemma  1 , 

0 -(n-2/  0 A 


but  since  [A  ,A]  = 0,  hence  [(A  )-,A]  = 0,  both  sides  of  the  equation  must  equal 

N N 

zero.  Here  QfH^)  just  depends  on  the  coefficients  of  and  their  der ivativ€*s . We 

N N 

thus  have  the  term  Q(HA)  is  an  exact  derivative  in  the  coefficients  of  H^,  i.e. 

N N N] 

Q(Hft)  = DQ' (H^) , and  in  fact  Q'(*)  is  easily  computable.  We  thus  have  = -Q' (H^) , 

from  which  we  deduce,  using  (4.1), 

(4-2)  Vn  = PnA  + *(HaJ'  *‘V  = <Vh  + JT  1(Po,n(tHA'A1?))  ' 

and  we  note  the  coefficients  of  ip  e h are  polynomial  functions  of  2H_  n<  i<  n-2, 

l 

and  their  derivatives.  We  view  ip  as  an  operator,  \p  : A . . . , -*  A . 

- (n-1) , “1  u,n-i 

We  now  bracket  (4.2)  with  A,  yielding 

(4.3)  [A,PN+nJ  = [A,PNJ  ' A + • 

dhn 

By  Theorem  4,  [A, Pi  = A = H (J - — ) , with  J = A - A the  matrix  differential  operator 

N Da 

which  determines  Hamilton's  equation  (see  Section  3),  and  H defined  in  the  obvious 


way,  i.e.  H]  I = A - £ . Also  note  that  the  right  hand  side  of  (4.3)  depends  only 


on  H^,  by  Theorem  4,  and  is  contained  in  n 2 s^-nce  the  left  hand  side  is.  We 

can  thus  define  the  (n  - 1)  x (n  - 1)  matrix  differential  operator  M by 

N N 

(4.4)  H W 52_)  = H{J  £2-)  • A + [A,0>  (h“)  ] , 

Da  Da  A 

DHN+n 

and  since  [A,P„  ) = J (- ),  (4.3),  (4.4)  implies 

N+n  Da 

DHN  DHN+n 

(4.5)  — > < 

Da  Da 

the  standard  form  of  Lenard  recursion  relations,  since  M,  J are  matrix  differential 


operators.  Equation  (4.2)  is  equivalent  to 


(4.6)  p = l * At_V,  <p  = *(h'V~  ’‘n+S),  0 < s < n - 1 . 

s+nt  v,s  v,s  A — — 

v=0 

We  observe  that  (4.4),  (4.5),  (4.6)  are  enough  to  establish  [A,P  ) = H(J HN)  , 

N A 

providing  we  define  P = 0,  N < 0,  and  we  shall  think  of  these  as  the  Lenard  relations. 
N 
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One  computes  M(j^)  = A(HAA)+  - (Allfl)+A  + ± ± (— ) A • (-iD)V‘1(PQ  ([H  ,A]C)). 

v=l  ‘ ' 

Also  note  equation  (4.5)  immediately  implies  Corollary  1,  Section  3. 

Remark  7 ■ The  self-adjoint  case,  and  for  that  matter  the  case  of  Remark  6,  may 

be  done  in  precisely  the  same  way  as  above,  but  has  some  novel  features.  For  if 
n-2 

A = £n  + £ (a.E3  + (£  - iD3a.),  a.  real,  then  if  <l>  is  a formal  eigenvector  with 

j=0  3 3 3 

eigenvalue  A of  the  associated  differential  operator,  we  have  by  the  self-adjointness 
of  A, 

(4.7)  • ~~  = (((iD)  3il)iJ/*+  ifr(  (iD)  3i|<)  *],  j = 0,1, ...,n  - 2 , 

dj 

where  * denotes  complex  conjugation.  And  so  (4.5)  implies 

(4.8,  «(£)  = ^ 

* 

Note  that  since  (M  > e I,  it  plays  the  rate  of  a constant  in  (4.8),  (see  beginning  cf 
Section  3) . 

Remark  8 . One  conjectures  from  computational  evidence  that  M in  fact  defines 
a symplectic  structure,  like  its  counterpart  J.  One  also  conjectures  that  (4.4), 

is  in  fact  a functional  identity  in  its  arguments,  i.e.  (4.4)  doesn't  depend  for  its 
. . N 

validity  upon  the  substitution  of  into  its  argument.  In  addition,  equations  (4.4), 

(4.5)  have  a bias,  at  least  in  their  derivations,  for  instance  if  instead  of  using 
N+l  N N N 

A = A • A,  we  had  employed  An  = A • An,  we  would  have  certainly  gotten  a 

different  tp,  but  we  doubt  that  we  would  have  gotten  a different  M;  however,  we 
cannot  prove  this,  and  merely  conjecture  this  to  be  so. 

2 

Examples.  For  the  self-adjoint  case,  A = -D  + aQ,  we  have 

J = D,  M = - i (D3  - 2 (a  D + Da  •) },  <Mv)  = i <2vD  - (Dv))  . 

4 0 0 4 

These  formulas  are  certainly  well  known.  For 

A = iD3  + i (a^D  + Da^ • ) + aQ  , 


— = < n> 

Da 


Tv  DX_ 
Da 


we  have 


r 


M = 


0,  D 
D,  0 


V M2 


■V  V 


M « j ID3  + axD  + Da^ • ) 


«2  = V + ! <%>  ' M3  " V + \ <Da0> 


M,  = 5-  t°5  + 5a, D3  + D3 (5a  ) • + (8a?  - 3(D2a,))D  + D(8a2  - 3(D2a,))  ■ ] 
49  1 1X1  11 


L 2 


= j {2iv2D2  + (Vj  - i (Dv^)  ) D + j (Ba^v^  + (D2V2  - Dv^  } 


These  formulas  for  J , M have  been  found  independently  by  H.  McKean  [6].  For 


4 2 2 

A = D + (a2D  + D a2‘)  + i(a^D  + Da^-)  + aQ  / 


then 


| 

Jn' 

J12' 

~-s 

J13 

Mll' 

M12' 

— ■ 

M13 

J = 

J21 ' 

J22' 

J23 

, M = 

M21‘ 

M22' 

M23 

, 

where 

J31' 

J32' 

J33 

M31' 

M32' 

M33 

J 

11  " J 

►o 

II 

(N 

21  = J23 

J32 

= 0 , 

J13  ■ J31  = D'  J22  " - I D'  J33  = - I ^ + (a2°  + DV>  ' 

2 

and  letting  D tp  = ip" , etc.,  we  have 

Mn  ■ ! °3  + I (a2°  + Dv>'  M12  = V21  ■ f aiD  + I ai  ' 

«13  = -^l  ■ -?  D'  + ? a2D3  + f ^ + (?a2  + VD+  <!a2"  +!ai>  ' 

M22  ■ I (|°5  + I [a2°3  + °\]  ‘ I ((a2  + a0  - a2,D  + °(a2  + "o  " a2>-11 
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M2  3 = '*32 


'V'  * I dlp2  ' (a'i  + 2 a2al’D  ' (2  aIa2  + 4 a]'"1  • 


33  8 


“ D7  + \ [(al  - a" 


1 a ) + D3(a^  - a"  - }-  a. ) • ] 


2 0 


2 2 0 


, , 1 i i . ■ ’2  32  1 , , , '2  32 

f I <t  a'  ' -a  + — a.  + a.a  )D  + D(-  al'  - a.  + - a + a_,a„) ) 
4 2 2 41  20  42  2 41  20 


wh  i 1 e 


i^(v)  = v 


2 V3p3  ~ 4 (V3  + iv2)c>2  + (vx  + iv2*  + ^iV31D 


\ t|  <-3vj  - iv“  + v'">  - 2$^'  + 3i^2v3  - i<‘1v2)  . 


Clearly  things  rapidly  get  out  of  hand. 

We  note  that  these  formulas  are  computed  from  (4.1)  — (4.4) . In  the  case 
n-2 

nr  1 

A = £ + l a^£  , one  can  use  the  above  formula  to  get  general  formulas  for  M,  J,  ip. 


i=0 


but  since  in  the  self-adjoint  case,  one  must  compute  the  dual  of  (^  - iD) J + = n_. 

by  hand  (up  to  order  5 n S , one  computes  M,  J,  ip  separately  for  each  n.  In 

[8],  Symes  computes  M,  J for  the  case  A = + a^5  + a^. 

Remark  9.  The  relations  of  this  Section  also  apply  to  the  systems  of  Section  2. 

We  derive  relation  (4.2)  for  these  systems  which  imply  relations  of  the  form  (4.4),  (4.5), 
(4.6)  . 

We  shall  use  the  notation  of  Theorem  1.  Observe 

N N+  N-  NO  N- 

L = (L  ) - (L  ) + (L  ) + 2<L  ) = PM  + B , 

N N 


where 


N + N - 
(L  ) - (L  ) , 


b = (ln)°  + 2(ln)~  = (h  ) + e , 

N N A m 


by  (2.14),  where  H = < — — — — AN  ^>,  and  6 e 4 . From  this  we  conclude 

N + 1 m -n,-m-l 
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r 


N + l 

p = L - B * (P  + B ) • L - B„  , 
N + l N+l  N N N+l 


= V + % • U - Bn+1 


‘V  + (Va 


L1  + 10m  • L - bn+i> 


hence,  remembering  (P  ,)  f A,  , we  have  shown 
^ N+l  1 , n 


(4.9) 


(pn+i>  ‘ (pnl)  + i(Va  ‘ L1 


(PL)  - A((H > > . 
N N A 


The  above  defines  the  function  A ( - ) . Taking  the  negative  transpose  of  the  above 

- - t T T 

equation,  we  have  P . = (LP  ) + A ( (H  ) ),  where  we  have  used  P„  = -P„,  I.  = L,  but 

N+l  N N A N N 

(4.10)  (P  L)  - (LP  ) = [P  , L]  = ( [A , ] ) T,  (by  Theorem  1)  , 

N -n  N -n  N -n  A + 

and  so  (4.9),  (4.10),  and  the  expression  for  P , , imply 

N+l 

p.ul  = PM  ' L + tM(H  > ) , 

N+l  N N A 

with  i|i(H  ) = (HftL)+  - (LH^)  + ((Ha,A)+)T.  Thus  a relation  of  the  form  (4.2) 

has  been  established,  as  claimed.  We  also  note  that  the  conjectures  of  Remark  8 are 

relevant  in  this  case. 


We  give  the  results  for  the  Toda  system,  i.e.  A. 


A.  .'s  equal  to  zero.  Define  E..=l  if  j = i + 1, 

10  , 13 

T / 3 

for  E.  ,.  Let  M,  J act  on  in  general  the  vector 
43 


H = 
a 


3H  3h 

3a.  3a  . 

0 n-1 


etc . for  H.  . Then  we  have 
b 


=b.,  A.  . , =a.,  all  other 
l i,i+l  l 

zero  otherwise  and  similarly 
, H = H (a,b)  , 


J 


0,0,,  the  n x n,  and  (n+l)  x (n  + 1)  zero  matrices  respectively,  and 
n n+l 

[S ] . . = (-6..  + E..)a.,  i = 0, . . . ,n  - 1,  j = 0, . . . ,n. 
l]  13  11  1 
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M = 


where 


AV  W3 


-MT,  M 

3'  2 


(M,  1 

a . a . m 

1 / _ T . • • 

= 1 (E..-E..),  1,3  = 0, 

1 

ID 

2 13  13 

[Ml 

2 T 

= 2a. (E..  - E..),  1,3  = 0,...,n 

2 

13 

1 13  13 

(Ml 

= a . (b.  E.  . - b.S..),  i = 0,.. 

3 

ID 

1 l+l  i3  1 13 

Finally  if  w = (uQ, 

— 

T 

.,u  ,,v,...,v)  , then  we  have 

n-1  0 n 

till  (w)  ] . . = \ & . . (a.u.  - a.  .u.  ) - a.  (E.  .v.  -ET.v  ),  i,  j = 0,1, . . 

i]  2 i]  l l l-X  l-l  l 13  l i]  l+l 
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